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Abstract 

Amphiphilic block copolymer solutions form various micellar structures including micelles 
and vesicles. We applied the density functional theory for block copolymers which we have 
proposed to amphiphilic block copolymer systems. The 3 dimensional simulation for AB 
diblock copolymer solutions and AB diblock copolymer / A homopolymer blends has been 
done and it is shown that the spherical micelles, cylindrical micelles and spherical vesicles are 
formed. It is also shown that the phase diagram for AB diblock copolymer / A homopolymer 
blends qualitatively agrees with the phase diagram obtained by the experiment. 

1 Introduction 

Block copolymers dissolved in solvent or other polymers assemble spontaneously to form various self 
organized structures such as spheres, rods and vesicles pQ. What structure is formed is a question 
very important in colloid science, but difficult to answer since the structure depends on many 
parameters such as the block ratio, the affinity between component units, the polymerization index 
and the volume fraction of the block copolymers. Attempt has been made to predict the micellar 
structure by computer simulation. Since the process of self-organization involves many molecules 
and is very slow, atomistic molecular dynamics simulation cannot be used for this purpose. 

Various coarse grained models have been used to study the formation of the micellar structure. 
Larson et al. 0EIE] used a lattice model to show that amphiphilic molecules actually form micellar 
structure . Bernardes |Hj studied the vesicle formation by the Larson type lattice model. Yamamoto 
[7] demonstrated that the DPD (dissipative particle dynamics) can show the formation and 
various structural change of vesicles. Similar attempt has been done by Noguchi El et al. for 
small molecule surfactant solutions. They studied formation and deformation of surfactant vesicles 
by Brownian Dynamics (BD). 

In these models it is not straightforward to connect the coarse grained model with the underlying 
atomistic model. (Some parameters can be associated with the mean field parameter. For example, 
the interaction parameters used in DPD can be mapped onto the Flory-Huggins \ parameter [rjlllUj. 
But this is not always successful, and needs more validation.) 

Other approach which takes into account of the molecular characteristics is the dynamic mean 
field theory which is based on the self consistent field (SCF) theory jTJElEUE]- In this approach, 
it is possible to determine the model parameters from the atomistic models for certain classes of 
polymers. 
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For periodic systems, computationally efficient and fast SCF method has been developed 15 
which uses eigen functions. The method has been applied for AB diblock copolymer / A homopoly- 
mer blends which forms micelles |16j . However the method is not best suited to study the micelles 
or vesicles which does not form periodic structures. 

For non-periodic systems, the real-space SCF simulation method has been frequently used. 
The method is based on the real-space calculation of the path integral and needs large CPU power 
and memory. Although considerable development has been achieved due to the improvement of 
algorithms an d in models it is still computationally demanding for large systems. 

Some works has been done for micellar systems using the real-space SCF simulation. For 
example, Fraaije et al. [T7] and Lam et al. ^2] have studied rather dense block copolymer solution 
with the mean field parameters calculated from the experimental data and shown that they form 
the micellar structure. Fraiije and Sevink |19j also demonstrated various intriguing self organized 
structure formed in the droplet of concentrated solution of block copolymers in equilibrium with 
the outer solvent phase. 

While these applications of SCF have achieved successful results, the dilute micellar systems 
have not been studied yet except for a few works. He et al. [201 - performed 2 dimensional SCF cal- 
culation for the dilute AB diblock copolymer solution (AB diblock copolymer / S solvent mixtures) 
and showed that the AB diblock copolymer actually forms 2 dimensional vesicles. Though this is 
an important achievement, the 2 dimensional simulation poses problems in the interpretation of the 
results. The spherical vesicles are 3 dimensional objects and cannot be handled by 2 dimensional 
simulation. At present, the 3 dimensional real-space SCF simulation for micellar systems costs too 
much CPU power and memory and difficult to perform. 

The computational cost of the dynamic mean field theory can be reduced drastically if one knows 
the free energy of the system as the functional of the density distribution of each component. The 
crucial point of this approach is the expression of the density functional. Earlier works (211122) used 
heuristic arguments for the expression, and was difficult to generalize. In the previous paper |23], we 
have proposed a general expression for the density functional which involves the same parameters 
appearing in the SCF theory. We have shown some preliminary results which shows the formation 
of the self organized structure for the mixture of block copolymers and homopolymers. 

In this paper, we shall show that our free energy functional can be applied and is effective for 
the study of 3 dimensional micellar systems. We shall show that micelles and vesicles are formed 
from the homogeneous state. (As far as we know, this is the first 3 dimensional simulation of the 
vesicle formation from the homogeneous state using the continuous field model.) We shall discuss 
the relationship between the micellar structure and the molecular parameters (volume fractions, x 
parameters, the block ratio). We shall also discuss the comparison between the present simulation 
and experiments and the SCF simulation by He et al. |2"H] . 

2 Free Energy Functional Model 

In this paper, we discuss two systems: one is the AB diblock copolymer solution (the solvent being 
a small molecule of C), and the other is a polymer blend, a mixture of AB diblock copolymer and 
A homopolymer. In both systems the block copolymers act as the surfactant and form the self 
organizing, micellar structures. 

In the density functional theory, the free energy F of the system is written as a functional of the 
density distribution of the monomers belonging to various subchains in the system. The subchains 
are indexed by (p, i), where p and i refer to the polymer, and the subchain in the polymer. In the 
present problem, the system consists of three subchains (AB,A), (AB,B) and (S,S), each of which 
stands for the A block of the AB diblock copolymer, the B block of the AB diblock copolymer, 
and the matrix, which is the C solvent in the case of solution, and A homopolymer in the case of 
polymer blend. 

Let 4>pi(v) be the volume fraction of the monomers belonging to the (p,i) subchain at point r. 
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We use the following free energy functional model |23) . 

F[{<t>pi{r)}} =Xj/ drdr' 2^/f~f~A Ptij g(r - r')^ pi (r)^ M -(r') 

+ E / dr fpi C p 

(r)hi(j) pi (r) 
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(1) 
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where Q(r — r') is the Green function which satisfies — V 2 5('" — r') = <5(r — r'), / P j is the block ratio 
of the i-th subchain of p-th polymer, b is the effective bond length and Xpi,qj i s the Flory-Huggins 
X parameter. A p ^j and C Pt ij are determined from 6, the degree of polymerization N p and the 
structure of block copolymer (see Ref. [231 f° r detail). 

The free energy functional eq Q reduces to the following form for the AB diblock copolymer 
/ S solvent blends. 
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Here we dropped the index p, which specify the polymer species, for simplicity. Thus (j> A {r), 4>b(t) 
and ips{r) stands for the density of the A component of the block copolymer, the density of the B 
component of the block copolymer and the density of the matrix (C solvent in the case of polymer 
solution and A homopolymer in the case of polymer blend). A^j and dj for diblock copolymer 
(i.j = A, B) is given in Appendix lAl 

It is convenient to introduce the new order parameter (the i/i-field) defined via the following 
equation [23; 

^(r) = v^W (3) 
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With this order parameter, eq is rewritten as 
F[ip A (r),ipB(r),ips(r)] 



E J drdr'2^Jj~A l] g(r~r')iP l (r)i; 3 (r') 

dr [2f A C AA ip A (r) lnV^(r) + 2] B C BB ^ B (r) ln^(r)] 
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3 Simulation 

We performed the simulation to get the equilibrium structure of the diblock copolymer solution. 
To get the equilibrium structure numerically, we minimized the free energy functional eq (@J by 
using the algorithm shown in this section. 

3.1 Constraints 

The minimization of the free energy was done under two constraints each representing the conser- 
vation condition and the incompressible condition. 



(5) 
(6) 



drVf(r) =V<t>i, (i = A,B,S) 

E #(»•) = 1 

i{= A ,B,S) 



where V is the volume of the system and fa is the spatial average of fa{r). fas is the average 
volume fraction of solvent and (j) A , cj) B are given by the following equations. 



t>A = f A faAB 
^B = fB^AB 



(7) 
(8) 



where cf> A B is the average volume fraction of the AB diblock copolymer. fa\s and fas satisfy the 
following equation. 

4>AB + = 1 (9) 

Constraints eqs (0, © give the following terms for the free energy functional 10} . 

Constraints [M^, feW, M?)} = E / dr \ ^ + «( r )l tyiM ~ ( 10 ) 



where Ai and k(t) are the Lagrangian multipliers which correspond to the conservation and in- 
compressible conditions, respectively. 
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3.2 Numerical Scheme 



The free energy was minimized by evolving the -0-field iteratively by using the steepest-descent 
method. 

4 n ™\r) = A(r) - to [ Mi (r) + \^{r) + «(r)Vi(r)] (11) 



where ui is positive constant and /j,i(r) is the chemical potential defined as 

5F[ip A (r),ip B (r),ips(r)} 



(12) 



By substituting eq into eq (|12|l . we get the explicit form of the chemical potential (See Appendix 
lUfor detail). 

We employed the time-splitting method and the ADI method |2U to update ip( new \r) from 
i/j(r) by eq ifTTl) . 
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^ new) (r) = ^ (4) (r)-c( 3 )A^ (4) 



(13) 
(14) 

(15) 

(16) 
(17) 



where w^"^ is positive constant and /4 is the chemical potential calculated from ip!- u '(r). The 
advantage of using the ADI method is that we can choose a large value for uj compared with one 
used in the explicit method which we employed in our previous work [23] . The evolution with 
eq Ijl7|l can be regarded as rescaling of the V'-field to satisfy the conservation law (this is analogy 
to the wave function in quantum mechanics; we normalize the wave function to get probability 
density). Thus eq l|17|) was modified as follows. 



, (now) / \ 



dr 



(18) 



We also evolved n{r) as well. 



K (ncw) (r) = K{r) - uj {i) [tf{r) - 4>i 



(19) 



Eqs l(TBl) - lfT^|) . l(TH|) and ifT^I) were solved iteratively to get the steady state structure. All sim- 
ulations were started from the nearly homogeneous state with small Gaussian white noise and 
the periodic boundary condition was used. The evolution was conducted until the ip field do not 
change anymore. This structure is not necessarily the equilibrium structure: in most cases it is a 
mctastable structure. As our purpose is to search various metastable structures, we did not pursue 
to find the real minimum of the free energy. 

The numerical scheme described above produces the equilibrium (or metastable) structure 
efficiently For example, we can get 3 dimensional AB diblock copolymer vesicles (lattice points: 
128 x 128 x 128, detail parameters are shown in the following section) starting from the homogeneous 
state in about 26 hours on a 3.0GHz Pentium 4 PC. 

For the systems with small volume fraction of the block copolymer and/or small x parameter, 
spontaneous formation of the self-organized structure from the homogeneous state was not always 
successful. It is difficult to get the micellar structure by our scheme especially for these systems. 
Thus we used the following procedure. 



5 



We started the simulations from the nearly homogeneous state by using rather large \ parameter 
(for example xab = 1) and/or large volume fraction (for example §ab — 0.1). This gives a state 
which has a large density fluctuation after several hundreds iteration steps. We then switched 
Xab an d 4>ab to the target value and continued the simulation for several thousand iteration steps 
(about 2000 to 4000 in this work) to get the equilibrium state. Though arbitrary and artificial 
it may look, such procedure was needed in order to get the phase separated structure. In fact, 
without the procedure, we could not reproduce the phase separation in the region where the x 
parameter or the volume fraction are small. 

4 Result 

4.1 AB Diblock Copolymer Solutions (AB Diblock Copolymer / C Sol- 
vent) 

First we show the result of the simulation for the AB diblock copolymer solutions. We set the A 
monomer is solvent-philic and the B monomer is solvent-phobic (i.e. the x parameter between the 
A and C monomers is small or negative, and that between the B and C monomers is large). 

Figure ^ shows the result of the simulation. Here the surfactant is a medium sized block 
copolymer (Nab = 20) with short solvent-philic part (/a = 1/3) and long solvent-phobic part 
(Jb = 2/3). The other parameters are N c — 1,4>ab — 0.1,0c = 0.9, xab = 1>Xbc = 1-75, and 
Xac are taken to be 0.5 (a), (b) and -0.175 (c). The system size is 486 x 486 x 486 and the 
number of lattice points is 96 x 96 x 96. 

The figure demonstrates that the density functional method gives self organized structures 
spontaneously starting from a nearly homogeneous initial state. Various micellar structures - 
spherical micelles, cylindrical micelles, open bilayers (not closed membrane structures) and vesicles 
are observed. Notice that the micellar structure is not the same: there is a distribution of micellar 
size as well as the micellar structure. We believe that real micellar system will also have such 
distribution. 

The morphological change shown in Figure ^ can be qualitatively understood as follows. As 
Xac decreases, the affinity between the solvent and the solvent-philic part increases. Therefore 
the interfacial area between the micelles and solvent tends to increase. As a result, the micellar 
structure changes from spheres to cylinders and then to vesicles. 

Figure shows the 1 dimensional density profile of Figure ^c). It is observed that AB diblock 
copolymer is strongly localized at the micellar structure (the vesicle) and the density of unimers 
are quite small. 

To investigate the effect of the architecture of block copolymers, we fixed the parameter except 
for the block ratio and performed several simulations with different block ratio. The parameters 
are set to N A b = 20, N c = 1,4>ab = 0.1,0c = 0.9, sxab = l,Xsc = 1-75, xca = -0.175, system 
size: 326 x 326 x 326, lattice points: 64 x 64 x 64 and the simulations have been done for block 
ratio f B = 0,0.1,0.2,0.3,1/3,0.4,0.5,0.6,2/3,0.7,0.8,0.9,1. The result is shown in Figure El It 
is seen that if the block ratio of the solvent-phobic subchain is too short( i.e., if fg < 1/3), the 
block copolymer acts as the solvent-philic homopolymer. Thus the system is homogeneous and 
the micellar structure is not formed. As increasing the block ratio of the solvent-phobic subchian, 
spherical micelles and the cylindrical micelles arc formed (the relation between the morphology of 
micellar structure and the block ratio is studied by Ohta and Nonomura for the strong segregation 
limit |25|). As increasing the block ratio further, the vesicles are formed. In this case the block ratio 
of solvent-phobic subchain is large, i.e. the block copolymer is so-called crew-cut block copolymer. 
If one increase the block ratio further, the block copolymer acts as the solvent-phobic homopolymer 
and causes macro phase separation to form droplets. 

4.2 AB Diblock Copolymer / A Homopolymer Blends 

We now study the AB diblock copolymer / A homopolymer blends. Like the case of solution, the 
AB diblock copolymer form micelles in the matrix of A homopolymers. 
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From the technical view point, simulation of polymer blends is easier than the simulation 
of polymer solutions since strong segregation is easily achieved in polymer blends for small x 
parameters: large x parameter makes the interface sharper and causes numerical stability problems 
at the interface. 

Figure 0] shows an example of the simulation result for a large system. Here the system involves 
128 x 128 x 128 lattice points (the system size is 646 x 646 x 646), and the parameters are set to 
be Nab = 20, Na = 10, /a = 1/3, fg = 2/3, xab — 1,<^ab = 0.1,0a — 0.9. In this case, vesicles 
can be observed more clearly than Figure [Jc). This is considered to be due to the change of the 
strength of segregation (the segregation of this system is stronger than the system of Figure ^c)). 

Various self organized structures can be formed depending on the parameters such as the block 
ratio /a, /b, the volume fraction 0ab and the x parameter xab ■ To study the relation between 
these parameters and the micellar structure, we conducted a simulation of smaller system which 
involves 64 x 64 x 64 lattice points. The polymerization index for the AB diblock copolymer and 
the A homopolymer is set to Nab = 20, Na = 10, respectively. 

Morphological change can be observed by chaining the parameters. Figure shows the effect 
of the block ratio /a, /b- Other parameters are set to 0ab = 0.1,0a = 0.9, xab = 1. If the 
B block is small (Figure Eta), /b = 1/3), the diblock copolymer forms a spherical micelle which 
consists of inner core made of B block and the outer corona made of A block. As the fraction of 
the B block increases, the spherical micelles become unstable, and vesicles consisting of a bilayer 
of the block copolymers is formed. (Figure [^b) , fg = 0.5). The thickness of the B block in the 
bilayer increases with the increase of B block ratio. Finally, when fg becomes equal to unity, the 
macro-phase separation between A homopolymer and B homopolymer is observed. (Fig. Eld), 
/b = 1). This behavior is just the same as one of the block copolymer solutions. Figure El 
shows the effect of 0ab, the volume fraction of diblock copolymers. Other parameters are set 
t° Ia = 1/3, /b = 2/3, xab = 1- As the volume fraction of the diblock copolymer increases, 
the micellar structures change from the spherical micelles to the cylindrical micelles and vesicles. 
Figure shows the effect of the x parameter between the monomer A and monomer B. Other 
parameters are set to /a = 1/3, fg = 2/3, 4>ab — 0.1, 4>a — 0.9. As the x parameter increases, 
(i.e., as the antagonicity between the monomers A and B increases), the structures change from 
the spherical micelles to the cylindrical micelles and then to vesicles. This behavior is just like the 
case of the volume fraction change. 

These results are considered to be consistent with the experimental results of Eisenberg et 
al.. [20 1271 They studied the morphology of diblock copolymers in mixed solvent (water 
plus dioxane) for various block copolymer concentration and for various solvent composition, and 
observed the change of the micellar shape from sphere to cylinder and then to vesicles as the water 
content of the solvent increases. If one regard the change of the water content as the change of 
the x parameter between the solvent and the subchain of block copolymer which forms core, our 
results are qualitatively in agreement with their experiments. 

In order to demonstrate the similarity between the experimental results and our simulations, 
we conducted the simulation in the parameter space of 0ab and xab- Fig. El shows the result. 
The figure is in qualitatively agreement with Fig. 6 in Ref. Note that our simulations was 

not done for large systems and the size of the simulation cell may affect the final micellar structure 
(the finite size effect or the effect of the periodic boundary conditions). 

5 Discussion 

The micellar structures, including vesicles can be formed for the AB diblock copolymer solutions 
by using the free energy model eq (JJJ . The diblock copolymer density is dilute in our systems and 
the diblock copolymer strongly localized at the micellar structures (i.e. the system is the strong 
segregation). Note that such systems are difficult to treat by the previous density functional ap- 
proaches. Our free energy model can handle these systems qualitatively correct and the simulation 
with our model is much faster than the real-space SCF simulation. 

The dynamic vesicle formation process, however, cannot be discussed by our numerical scheme 
because our scheme does not satisfy the local conservation of mass. This is in contrast to the 
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particle methods (BD 8 and DPD 7.)- Nevertheless it is instructive to study the process of the 
micellar formation based on our model. 

The vesicle formation process for the AB diblock copolymer / A homopolymer blend by our 
scheme is shown in Figure El (N AB = 20, N A = 10, f A = 1/3, f B = 2/3, xab = 1, 4>ab = 0.1, 4>a = 
0.9, system size: 326 x 326 x 326, lattice points: 64 x 64 x 64, a/ 1 ) = 0.5). It is observed that 
the nucleus of micellar structure is first formed by the association of block copolymers (the matrix 
polymer is expelled from the micellar structure at this stage). The micellar structure then grows 
taking solvent into its core from the surrounding, and finally the vesicle is formed. 

This process is the same as that observed by He et al. [2H in their 2 dimensional simulation using 
the SCF method. (They also used the non-conserving scheme for the evolution of the local density.) 
According to the scenario of He et al., the nucleus of micellar structures grows first. If the nucleus 
grows sufficiently large, the nucleus takes solvent into the core because the hydrophilic subchains 
localized at the large micelles are energetically unfavorable. Finally the micellar structure changes 
their shape from the sphere-like structure into shell-like structure and thus vesicles are formed. 

Therefore our vesicle formation process is the same as one of He et al.. This agreement implies 
that the schemes which does not satisfy the local conservation can form vesicles easily compared 
with the schemes which is based on the dynamics and satisfy the local conservation (for example, 
the Cahn-Hilliard type TDGL equation or the dynamic SCF method J3j). The dynamic vesicle 
formation process based on the continuous field model simulation is the future work. 



6 Conclusion 

We have shown that the self-organized structure of block copolymers in solvent can be predicted 
by the density functional theory using the free energy proposed in Ref. 23 . This is the first 
3 dimensional simulation of vesicle formation by using the continuous field model. Since the 
theory can take into account of actual structure of polymers (degree of polymerization, block ratio, 
topological structure), it will be useful to understand and to predict the micellar structure of 
surfactant systems. 



Appendix 

A Form of Matrices Aij, dj 

The matrices Aij , Cy for diblock copolymer is represented as follows [231 ■ 



.4 



' N\ B Vf\P B 



B I . i 



SaSb f 2 



Nai 



^'>(\/-)-^vf 



4/A/i 



4/a/e 



(20) 
(21) 



where S 1 (£) is the inverse matrix for the normalized scattering function matrix S(£) for the ideal 
systems. 

J2Si j (0Si k 1 (0=Sik (22) 

3 



S(£) is defined as follows. 

4(0 = 



A(e-/-« - 1 + f A £) H (e-^ - l)(e-'*e - 1) 

1 e e 2 

-(e-^ - l)(e-'*« - 1) ^(e- fBi - 1 + M) 



(23) 
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B Chemical Potential 

The chemical potential for the system can be calculated by substituting eq Q into eq (|12|l . The 
chemical potential for the A, B subchains and S solvent is as follows. 



j(=A,B) 

+ 2f A C AA [2^ A {r)\n^ A [r)+ip A {r)] 



+ 4V ' fAfBC AB i>B{r) (24) 
+ J2 2 XAjtf{r)iJj A (r) 

j(=A,B,S) 



Mr)= J2 [ dr' 4^M~A Bj g(r - r')^(r') 

j(=A,B) 

+ VbCbb [2i>B (r) In ip B ) + tpB (r)} 



+ 4y r f A fBC AB ^ A (r) (25) 
- j V 2 ^ S (r) 

+ E 2 XB,jil>j(r)i>B{r) 

j(=A,B,S) 
2 

Ms(r) =^r- [2^>s(r) ln^g(r) + ips(r)] 
l\ s 

b 2 

- yV 2 feW (26) 
+ J2 2 XB, 3 tf{r)ip s (r) 

j(=A,B,S) 

Note that /Xj(r) contains the Laplacian term —b 2 V 2 t/ji(r)/3 for all i. Thus the evolution equation, 
eq (THl) . can be expressed as 

7 2 

# ew) (r) = &(r) + ^-V 2 Mr) + ■■■ (27) 

As mentioned before, eq H27[) can be solved stably by using implicit scheme (the ADI scheme, in 
this work), in analogy to the diffusion equation. 
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Figure 1: Micellar structures of the AB diblock copolymer solutions (AB diblock copolymer / C 
solvent blends) with different \ parameters. The system size is 486 x 486 x 486 (96 x 96 x 96 
lattice points). Parameters are set to Nab = 20, Nc = 1,4>ab = 0.1,0c = 0.9, xab — 1,Xbc = 
1.75, f ' A = 1/3, /s = 2/3. (a) xac = 0.5, (b) xac = 0, (c) xac = -0.175. The gray surfaces are 
isodensity surfaces for 0s(r) = 0.5. 
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Figure 2: Density profile of an AB diblock copolymer vesicle. The parameters used are the same 
as Figure H^c). The 1 dimensional density profile (right) is the profile along the black line in the 3 
dimensional isosurface data (left). 
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Figure 3: Phase diagram of AB diblock copolymer solutions, (a) fs — 1/3, (b) /s = 0.5, (c) Jb = 
0.6, (d) f B = 1. N AB = 20, N c = IJab = 0.1,0c = 0.9, X ab = 1,Xbc = 1-75, X ca = -0.175, 
system size: 326 x 326 x 326, lattice points: 64 x 64 x 64. 
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Figure 4: A result of the simulation for the AB diblock copolymer / A homopolymer blend. The 
system size is 646 x 646 x 646 which involves 128 x 128 x 128 lattice points. The parameters 
of polymers are Nab — 20, Na — 10, Ja — 1/3, /s = 2/3, xab = 1- The volume fractions are 
4>ab — 0.1, <j)A — 0.9. The gray surfaces are isodensity surfaces for ^s(r) = 0.5. 
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Figure 6: Results of simulations with different volume fractions. The system size is 326 x 326 x 326. 
N AB = 20, N A = 10, = 1/3, Sb = 2/3, xab = 1. (a) 4>ab = 0.025, (b) 4> AB = 0.0525, (c) 
4>ab = 0.075, (d) (j)AB = 1- The gray surfaces are isodensity surfaces for 0b (r) = 0.5. 



16 



Figure 7: Results of simulations with different \ parameters. The system size is 326 x 326 x 326. 
N AB = 20, N A = 10, Sa = 1/3, f B = 2/3,4>ab = 0.1, ^ = 0.9. (a) X ab = 0.0475, (b) X ab = 
0.525, (c) xab = 0.75, (d) xab — 1- The gray surfaces are isodensity surfaces for 0b (r) = 0.5. 
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Figure 8: Phase diagram of AB diblock copolymer / A homopolymer blends. The plotted symbols 
(circles, squares, etc.) correspond to the result of simulations with various xab and 4>ab- The 
bilayer phase means open, disk-like micelles (actually it is observed by experiment \27\). 



18 



5 10 15 20 25 30 




(C) 







1 

0.8 
0.6 
0.4 
0.2 


1 

0.8 
0.6 
0.4 
0.2 





10 



15 20 25 



30 



10 



15 

X 



20 



25 30 



Figure 9: Snapshot of vesicle formation process for the AB diblock copolymer / A homopolymer 
blend. Number of iteration step = 0, 470, 500, 530, 990 for (a), (b), (c), (d), (e), respectively. Nab = 
20,7V A = 10, /a = 1/3, /b = 2/3, \ab = h<t>AB = 0.1,j> A = 0.9, system size: 326 x 326 x 326, 
lattice points: 64 x 64 x 64 (The parameters are the same as Figure 0] except for the system size.) 
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